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ABSTRACT 

We present an equation of state and radiative opacities for a strongly magnetized hydrogen plasma 
at magnetic fields B, temperatures T, and densities p typical for atmospheres of isolated neutron 
stars. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and 
Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for 8x 10 11 
G < B < 3 x 10 13 G, 2 x 10 5 K < T < 10 7 K, and a wide range of p. We show that bound-bound and 
bound-free transitions give an important contribution to the opacities at T < (1 — 5) x 10 6 K in the 
considered range of B in the outer neutron-star atmosphere layers, which may substantially modify 
the X-ray spectrum of a typical magnetized neutron star. In addition, we re-evaluate opacities due to 
free- free transitions, taking into account the motion of both interacting particles, electron and proton, 
in a strong magnetic field. Compared to the previous neutron-star atmosphere models, the free-free 
absorption is strongly suppressed at photon frequencies below the proton cyclotron frequency. The 
latter result holds for any field strength, which prompts a revision of existing theoretical models of 
X-ray spectra of magnetar atmospheres. 

Subject headings: equation of state — magnetic fields — plasmas — stars: atmospheres — stars: neutron 



1. INTRODUCTION 

Models of neutron star atmospheres are needed for in- 
terpretation of their spectra and cooling. These atmo- 
spheres differ from the atmospheres of ordinary stars be- 
cause of the high gravity and magnetic fields (for review , 
see, e.g.. lPavlov et alJll995HVentura fc Potekhinll2001h . 

A magnetic field is called strong if the electron cy- 
clotron energy huj cc = heB/m c c exceeds 1 a.u. - i.e., 
the field strength B is higher than Bq = m1ce 3 /h 3 = 
2.3505 x 10 9 G, where m c is the electron mass, e the 
elementary charge, and c the speed of light. Usually 
the field is called superstrong if hjj cc > m c c 2 , that 
is B > B r = rale 3 /eh = 4.414 x 10 13 G. Most of 
the radio pulsars have magnetic fields B ~ 10 12 — 
10 13 G (jTavlor. Manchester, fc Lvnd 119931) . whereas 
anomalous X-ray pulsars and soft ga mma repeaters are 
thought to have supers trong fields (e.g. . lMereghettill2001t 
iThompson et al.l I2000L and references therein). Non- 
negligible amount of neutral atoms can exist in the pho- 
tos phere at typical neutron-star tempera tures T ~ 10 6 
K (jPotekhin, Chabrier. fc ShibanovH"l999l . hereafter Pa- 
per I). A strong magnetic field enhances atomic binding 
and makes the quantum-mechanical characteristics o f an 
atom dependent on its motion across the field (see lLa.il 
|2001| for a recent review). In photospheres of the neu- 
tron stars, the field is, as a rule, strongly quantizing, i.e., 
it sets all the electrons on the ground Landau level. This 
occurs if /3 e 3> 1 and p < ps, where 

& = hw cc /k B T » 134.3 B 12 /T 6 , (1) 

p is the density, and ps = m^/ (n 2 \/2 a 3 ^) w 7100 B^ 2 



3 (for the hydrogen plasma). Here and hereafter, 
r7i n + m n , m p is the proton mass, et m = (hc/eB) 1 / 2 



g cm 

777 H — "*p 

is the magnetic length, k B is the Boltzmann constant, 
B12 = B/10 12 G, and T 6 = T/10 6 K. 

Opacities for the two polarization modes of radia- 
tion are quite different in strongly magnetized plasmas 
(e.g.. iPavlov et all 119951 . and references therein), which 
makes ther mal emission of ne utron stars polarized and 
anisotropic (|Zavlin et al.lll995l ). T he mean opacities are 
stron gly reduced at e ^> 1 (e.g., ISilant'ev fc Yakovlevi 
1980); thus the bott om of the photospher e is shifted to 
high densities (e.g., lLai fc Salpeterl fl997t IPavlov et al.l 
[19951) . 

The chemical composition of neutron-star atmospheres 
is not precisely known. Just after the neutron star birth 
in a supernova explosion, the outer stellar envelope is 
most probably composed of iron. However, light ele- 
ments may be brought to the surface later (e.g., by fall- 
back, accretion, or encounters with comets). Because of 
rapid gravitational sediment ation, the lightest element 
will cover the surface (see iBrown. Bildsten. fc Chang 
120021) . About 10 12 — 10 14 grams of hydrogen (< 10" 19 
Mm) is sufficient t o fill t he entire photosphere. 

Shibanov et al. (1992) presented the first model of hy- 
drogen atmospheres with strong magnetic fields. Later 
it was developed beyon d the diffusion approximation 
([Shibanov fc Za vlin 199 3) a nd use d for astrophysical pre- 
dictions (e.g., Zavlin et a ll 119951; IHo fc Lail l200l 120031: 
Lai fc Hdl2^lOzelll200l[2^0l[ T Zane. Turolla. fc Trevid 



20001: I Zaire et al.l l2001l )and for interpretation of ob- 



served neutron-star spectra (e.g..lOze l. Psaltis, & Kaspi 
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fl995h . 

The above studies assume that the atmosphere is 
fully ioniz ed. Meanw hile, it was recognized long 
ago (e.g., iMillerl Il992f ) that a significant contribu- 
tion to the opacities of neutron-star photospheres 
with strong magnetic fields might come from bound- 
bound and bound-free absorption by atoms. Exam- 
ples of monochromatic opacities i n par tially ionized 
iron (iRaiagopal. Romani. fe Millerl 19971) and hydro- 
gen (|Potekhin. Chabrier. fe Shibanovl2000f ) atmospheres 
confirmed this conjecture. In Paper I we have presented 
an equation of state (EOS) of a partially ionized hydro- 
gen plasma for the values of T and B typical for at- 
mospheres of the radio pulsars. Here we report results of 
extensive calculations of thermodynamic functions based 
on the theory developed in Paper I, supplemented by 
calculations of the opacities (monochromatic and Rosse- 
land mean). Partial ionization and plasma nonideality 
are taken into account for 11.9 < log 10 B/G < 13.5 and 
5.3 < log 10 T/K < 7.0. Bound-bound and bound-free 
radiative transitions are treate d within the framework of 
a previously developed th eory (|Pavlov fe Potekhinlll995l : 
iPotekhin fc Pavlovlll997h . The free- free absorption cross 
sections are re-evaluated. Whereas the previous authors 
considered photoabsorption by an electron scattered off 
a fixed Coulomb center, we take into account the finite 
proton mass, which has a nontrivial effect on the pho- 
toabsorption in a quantizing magnetic field. 

The paper is composed as follows. In Sect.[2]we formu- 
late the main assumptions and give the basic formulae 
used in our work. Section presents the EOS of par- 
tially ionized hydrogen under conditions in neutron-star 
photospheres. In Sect. 0] we discuss various contributions 
to the hydrogen photoabsorption cross sections in strong 
magnetic fields and derive a new formula for the free-free 
cross section. Opacities of hydrogen photospheres of the 
neutron stars are discussed in Sect. [5] Appendices give 
some detail of calculation of the free-free cross sections. 



2. BASIC EQUATIONS AND PHYSICS INPUT 

2.1. Hydrogen Atom Moving in a Magnetic Field 

If an atom rests without motion in a strong mag- 
netic field, there are two distinct classes of its quan- 
tum states: at every value of the Landau quantum num- 
ber n and the magnetic quantum number —s (n > 0, 
s > — n), there is one tightly bound state, with binding 
energy growing asymptotically as [ln(_B/_Bo)] 2 , and an 
infinite series of hydrogenlike states with binding ener- 
gies approaching the ener gies of a field- free H atom (e.g., 
ICanuto fc Ventural Il977h . The atom is elongated: its 
size along the magnetic field B either decreases logarith- 
mically (for the tightly bound states) or remains nearly 
constant (for the hydrogenlike states), while the trans- 
verse radius is close to a m , decreasing as B -1 / 2 . The 
radiative transition rates are different for the three basic 
polarizations: the linear polarization along the field and 
the two circular polarizations in the transverse plane. 

This simplicity is destroyed when atomic motion is 
taken into account. The electric field, induced in the 
comoving frame of reference, breaks down the cylindri- 
cal symmetry. In the nonrelativistic quantum mechanics, 
the binding energies and wave functions of the H atom 
are given by a solution of the Schrodinger equation with 



the two-particle Hamiltonian 



H 



2m p 2m c 



(2) 



where mj, r^, and 7r^ are the mass, radius, and kinetic 
momentum of the electron [i =e) or proton (i =p) . The 
kinetic momentum (rela ted to the velocity) equals (e.g., 
lLandau fc Lifsiu^[l9"76h 



TTj = Pi - —A(ri), 

c 



(3) 



where qi is the charge of the ith particle (q c = —q p = 
— e), Pi is the canonical momentum (i.e., pi — —ihVi in 
the coordinate representation), and A(r) is the vector 
potential of the field. A conserved quantity related to 
the center-of-mass motion is the pseudomomentum 



K = 7T p + 7T C + -B x (r P 



(4) 



Let the z axis of the Cartesian coordinates (x,y,z) 
be dire cted along B. Separating the center-of-mass 
motion dGor'kov fc Dzvaloshinskiilll967t IPotekhin! 11994 
iVincke. Le Dourneuf. fc Bavel 119921 ) and choosing the 
gauge of the vector potential in the form 



A(r) = -Bx r ^ r 

2 V to h 



(5) 



where ro is arbitrary, one comes to the effective one- 
particle Schrodinger equation 



H ± + H K (r ) 



\ro 



ip(r) = Eip(r). (6) 



ro is a "shifted" relative coordi- 



Here, r 

nate, p z = —ihd/dz is the z component of its conjugate 
momentum p, 

^2 

(7) 



Hj_ = — —B ■ (r x p), 



2// 



m D c 



H K (r ) = -!—(K + -Bxr ) 2 
2mn V c 



K+-Bxr )-(Bxr), (8) 

mnc v c 



and \x = m e m p /mH is the reduced mass. In Eq. I|7|). 



it = p H B x r, 

2c 



0) 



and the subscript '_L' denotes a vector component per- 
pendicular to B. 
iJjc(ro) turns to zero, if we set ro = r c , where 

r c = x K (10) 

is the relative guiding center (the difference between the 
electron and proton guiding centers). This choice of ro 
is most useful for bound states with large K± and the 
states of the continuum, whereas for bound st ates with 
small K i the choice r n = is m ost appropriate (|Potekhinl 
fl99l IPotekhin fc Pavlovlfi997l) . 

The eigenfunctions of H± are th e Landau functions 
&ns(r ±) (given, e.g., by Eq. (5) of IPotekhin fc Pavlovi 
(1993), with cigcncncrgics 



E, 



nhuo cc + (n + s)huj C p, 



(11) 
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where w cp — (m e /m p ) Lo ce is the proton cyclotron fre- 
quency. 

It is convenient to expand the wave function in the 
basis of ® ns (r±), 



(12) 



and to label tp{ r ) by numbers n and s, corresponding to 
the leading term of this expansion. The third quantum 
number v then enumerates "longitudinal" energy lev- 
els. The adiabatic approximat i on widely used in the past 
(e.g^lCanuto fe Venturalfl977l ; lGor'kov fe Dzvaloshinskiil 
I1967D corresponds to retaining only one term n! = n, 
s' = s in Eq. (H2). We perform calculations without this 
approximation. 
The total energy of the atom in Eq. ([6]) can be written 

as 

E = E±+El v (K ± ). (13) 

Here, the longitudinal energy e\ sv {K_l) is negative for 
the bound and autoionizing (resonance) states and posi- 
tive for the continuum states, in which the motion along 
z is infinite (in the latter case, v is continuous). Since 
n = for the bound states of H atom in a strong mag- 
netic field, we will drop the number n but imply n = 
for these states. Then the binding energy is 

e,„(Ax) = \El(K x )\ - shiu cp . (14) 

The substitution of Eq. (fi"2"j) in Eq. © reduces the 
problem to the set of the coupled channel equations 



(p 2 z /2ti + E± s ,-E)g n , s ,(z) 



where 

is a total coupling potential, and 



V n . 



'(ro,z) 



V n . 



<(r ,z) = (ns\ - e 2 /\r Q 



'I \n's') 



(15) 

(16) 
(17) 



is an effective Coulomb potential. Here, 
(ns\f(r)\n's') ± = /"*™(r±)/(r)*„^(rx)d 2 r x . (18) 



Numerical s olutions o f Eq. J6l) for various K± were 
presented by IVincke et al.1 (|1992f L At superstron) 
fields, binding energies were calculated bv lLai fe Salpete 
(1995). The system of equations (fl~5|) was numericall y 



solved for the discrete atomi c states by iPotekhinl (|1994j) . 
and for the continuum by Potckhi n fc Pavlovl l|1997[ ). 
According to these studies, an atom moving across the 
strong magnetic field acquires a constant dipole moment 
parallel to r c . Those radiative transitions, which were 
dipole-forbidden for an atom at rest because of conser- 
vation of the z-projection of the angular momentum, be- 
come allowed and should be taken into account in the 
atmosphere models. If K± is small enough, the dipole 
moment is also small. When K± exceeds a certain criti- 
cal value, the atom becomes decentered: the average dis- 
tance between the electron and proton approaches r c . In 
this case, K± characterizes the electron-proton distance, 
rather than the atomic velocity. The binding energies 



(fT4|) decrease with increasing K±. Asymptotically, at 
large K±, all longitudinal energies tend to — e 2 /r c . In 
this limit, the cylindrical symmetry of the wave function 
and dipole selection rules are restored, but the axis of 
symmetry is shifted to the distance r c from the Coulomb 
center. 

2.2. Thermodynamic Model 

The EOS for partially ionized hydrogen in strong mag- 
netic fields was constructed and discussed in Paper I. 
We employ the free energy minimization technique in 
the "chemical picture" of a plasma ( for discussion of 
its advantages and limitations see, e.g ., |Potekhinlll996bl : 
iSaumon. Chabrier. fc Van Hornlll995f) . The treatment is 
based o n the framework of the free energ y model devel- 
ope d bv ISaumon fc Chabrierl (ll99lL|1992f) at B = (see 
also lSaumon. Chabrier. fc Van Hornll 19951 and Sect. II of 
Paper I) and extends it to the strong magnetic field case. 

We consider a plasma composed of N p protons, N c 
electrons, Ah hydrogen atoms, and A mo i molecules in a 
volume V, the number densities being nj = Nj/V. The 
Hclmholtz free energy is written as the sum 



F = Ft 



F n 



(19) 

where F? d , F? d , and F^ n are the free energies of ideal 
gases of the electrons, protons, and neutral species, re- 
spectively, takes into account the Coulomb plasma 
nonidcality, and F™ u is the nonideal contribution which 
arises from interactions of bound species with each other 
and with the electrons and protons. In Eq. (fl9|) we have 
disregarded the additive contribution due to photons, 
since it does not affect ionization equilibrium. Moreover, 
generally we need not to assume thermodynamic equi- 
librium of radiation with matter. Ionization equilibrium 
is given by minimization of F with respect to particle 
numbers under the stoichiometric constraints, provided 
the total number Nq of protons (free and bound) is fixed. 
The latter number is determined by the total mass den- 
sity: n = Nq/V w p/tor. 

The first term in Eq. HSJ is Ff d = (J, e N e - P C V, where 
jU e and P e are the chemical potential and pressure of the 
ideal Fermi gas, respectively. They are obtained as func- 
tions of the electron nu mber density n e and temperat ure 
T from equations (e.g.. iBlandford fc Hernquist]ll982f ) 



br, 



3/2 oo 



= P T -J-r V(2 - 5 n0 )(l + 2bn) 1 /% /2 ( Xn ,r n ), (20) 

v 2tH fr^ 

2~ d 



n=0 
oo 



n c = A 3 



n=0 



T n dXn 



h/2(Xn,T n ) 



,(21) 



where 



h/ 2 (x,T)= (e t -x + l)- 1 ^t(l+rt/2)dt, 
Jo 

T/T x 

Xn — u rr, + T 



-1 _ -1 



y/l + 2bn ' 



and b = (B/B r ). In these equations, Ac = h/(m c c) 2 as 
3.8616 x 1CT 11 cm, P r = m c c 2 /A c w 1.4218 x 10 25 
dyn cm -2 , and T r — m e c 2 /k& ~ 5.930 x 10 9 K are 
the relativistic units of length, pressure, and tempera- 
ture, respectively, and m e c 2 is not included in jjb c . We 
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employ analytic fitting formulae to the standard Fermi- 
Dirac integral I\/2 in Eq. (f20|) and to the expression in 
square brackets in Eq. (|21|) . accurate within a few parts 
in 10 3 . These fits are pre s ented , respectively, in Sect. 
Ill of IChabrier fc Potekhinl (|1998l ) and in Appendix C of 



iPotekhinl (|1996aT L When the electrons are nonrelativis 
tic (as usually in the photospheres), Eqs. (|20|) and ([2Tjl 
reproduce Eqs. (30) and (31) of Paper I. 

The Coulomb free energy contribution consists of three 
parts, F^ = F pp + F cc + F pe , which represent, respec- 
tively, the proton-proton, electron-electron, and proton- 
electron interactions. There was no detailed study of the 
influence of a strong magnetic field on these contributions 
in the p and T domain we are inter ested in. Therefore we 
employ nonmagnetic expression s (|Chabrier fc Potekhinl 
ll998tlPotekhin fc Chabrierll2000D . scaled with B. Specif- 
ically, the nonmagnetic expression for F pe is based on 
numerical results obtained in the hypernetted chain ap- 
proximation for the linear resp onse theory with a local 
field correction (|Chabrierl [19901 . The nonmagnetic ex- 
pressions for F pp and F ee result from fitting the most 
accurate numerical results available in the literature (see 
iPotekhin fc Ch abricr 2000 for references). In the strong 
magnetic field, the -B-scaling of the nonmagnetic F^. is 
constructed so as to match known low- and high-density 
limits (Sect. IIIB of Paper I). 

The free energy of nondegenerate and nonrelativistic 
gas of protons, F id , is given by 



FP d /N p k B T = ln(2TTal L X p n p ) + In (l - e^) - 1 

+/3 p /2-ln[2cosh( 3p /3 p /4)], (22) 

where g p = 5.585 is the proton gyromagnetic factor, A p = 
H \/2ti I (k B Tm p ) is the proton thermal wavelength, and 
/3 P = hLu cp /k B T « 0.0732 B 12 /T 6 . 

Let N sl/ be the total number of H atoms with given 
quantum numbers s and v in the volume V, and let 
p S v(K±) d 2 K± be the probability for such atom to have a 
transverse pseudomomentum in an element d 2 K± around 
K± . Then the ideal part of the free energy for hydrogen 
atoms is 



F id - = k b ± 



In 



(2^) 2 
/Ah — y — p s 



■1 - e su (K x )/(k B T) } p„(K ± ) d 2 K ± 



+N B k B T {(3 p /2 - In [2 cosh( 5p /3 p /4)] } , (23) 

where Ah ~ A p is the thermal wavelength of an atom. 
The probability density p sv {K±) is calculated in a ther- 
modynamically consistent way from derivatives of the 
total free energy F with respect to the particle num- 
bers. Molecules H2 are treated in an approximate man- 
ner, without taking into account their excited states and 
possible effects caused by their motion across the mag- 
netic field and rotation. Finally, the nonideal part of 
the free energy of neutral species, F^° u , is obtained in 
frames of the hard-sphere model, with effective radii de- 
pending on the quantum numbers and pseudomomenta 
of interacting atoms (see Paper I for detail). 

Once the free energy is obtained, its derivatives over p 
and T and their combinations provide the other thermo- 



dynamic functions. 



2.3. Polarization Modes and Opacities: Basic 
Relations 

Propagation of radiation in magnetized plasmas 
was discussed in many papers and monographs (e.g., 
lGinzburg|[l970h . At photon energies hw much higher 
than 



pi 



47rfi 2 e 2 n e 



m c 



1/2 



28.7pJ /2 eV, (24) 



where uj p \ is the electron plasma frequency and po = 
p/g cm -3 , radiation propagates in the form of extraor- 
dinary (hereafter labeled by index j • = 1) and ordi- 
nary (j = 2) normal modes. These modes have dif- 
ferent polarization vectors e-,- and different absorption 
and scattering coefficients, which depend on the an- 
gle 9b between the propagation dire ction and B (e.g., 
iKaminker. Pavlov, fc Shibanovl Il982f ) . The two modes 
interact with each other via scattering. Vec tors e, for 
a fully ionized plasma have been derived by IShafranovl 
(|1967f h IVentural (|1979f ) gave an instructive analysis of the 
plasma polarization mod es relevant to the neutron stars. 
iGnedin fc Pavlovl |l973) formulated the radiative trans- 
fer problem in terms of these modes. They introduced 
the convenient real parameters q and p, which completely 
determine the normal mode polarization properties, and 
which are defined as 



q + ip 



e xx cos 2 6 B + e xz sin 28 B - s zz sin 2 9 B 
2 i (e xy cos 8 B + e yz sin 8 B ) 



where E p are the compo nents of the complex permittivity 
tensor (|Ginzburg|ll970f L and the z axis is directed along 
B. The parameter q determines the ellipticity of the 
normal modes, and the parameter p is associated with 
absorption of radiation. In the most common case, one 
has 

M»|p|, g( w ,0 B )«g( W )£^-. (26) 

Z COS o B 

These relations may be invalid in narrow frequency 
ranges where resonant absorption occurs (e.g., near the 
electron or ion cyclotron resonance). 

The formulae for ej which take in to account contri- 
bution of the plasma ions, impli ed in IShafranovl (fl967). 
have been explicitly written by iHo fc Lail (|20'oil ). The 
electron-positron vacuum polarization in a strong mag- 
netic field dramatically changes the n ormal-mode prop- 
erties in certain p-uj domains (see iPavlov fc GnedirJ 
11984 for a review) . The vacuu m dielectric tensor has 
been obtained by Adler ( 19711 ) at B <C B T and by 
iHevl fc Hernquistl (tl9J7jEl) in bo th limits o f B < B T 
and B » B r . Using these results, IHo fc Lail (|2003f ) de- 
rived convenient expressions for the polarization vectors 
of normal modes, which take into account the contribu- 
tions from the electrons, ions, and vacuum. 

The presence of bound species modifies the com- 
plex permittivity tensor and hence the properties of 
the normal modes. Their accurate treatment in a 
partially ionized medium with a strong magnetic field 
is a complicated problem, which has not been solved 
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yet. The normal polarization modes of a neutral 
gas of hydr ogen atoms in s t rong m agnetic fields were 
studied by iBulik fc Pavlov! (|1996l ). who applied the 
Kramers-Kronig relations to the bound-free and bound- 
bo und atomic absorptio n coeffi cients obtained previously 
bv iPotekhin fc Pavlovl (|1993h and iPavlov fc Potekhinl 
(|1995h . The thermal motion effects, which had not been 
calculated by that time for the bound-free transitions in 
strong magnetic fields, were evaluated using a pertur- 
bation approximation (|Pavlov fc Meszarosl [l993). The 
qualitative behavior of the polarization vectors proved to 
be the same as for the fully ionized plasma in a wide range 
of uj and Ob, where Eq. (|26|) holds. However, there are 
quantitative differences. The frequencies at which Eq. 
(|26|) is not valid, are of the order of the photoionization 
threshold lu^ and the principal bound-bound transition 
frequencies, where the resonant absorption takes place in 
the neutral gas (instead of uj cc and w cp in the fully ion- 
ized plasma). The numerical values of q(u>) are modified. 
For instance, at ui well above ujth, <z( w ) is larger than it 
would be in the case of full ionization, which makes the 
normal mode polarization more linear. 

Polarization properties of normal modes in a partially 
ionized plasma with a strongly quantizing magnetic field 
remain unexplored. Hereafter we consider mostly situa- 
tions where the neutral fraction is small. Since even for 
completely nonionized gas the properties of the polar- 
ization vectors (in particular, their dependence on 9 b) 
are qualitatively the same as for the fully ionized plasma 
(except for the resonant absorption frequency ranges), we 
assume that the polarization modes for the fully ionized 
plasma a re a good appro ximation and adopt the formulae 
given bv lHo fc Lail pOOl . 

At a fixed photon frequency uj, the absorption opacity 
Kj(8b) in each mode j and scattering opacities k s ^,{9b) 
from mode j into m ode j' can be presented as (e.g., 
iKaminker et al.lll982D 

l 

K^e B )=m^ |ej,a(0n)| 2 oS. (27) 

a=-l 
a=-l 

x— f \e f , a (e' B )\ 2 sind' B dd' B , (28) 

"»H Jo 

where a = 0, ±1, e^o = e JiZ is the ^-component of ej, 
andej ; ±i = (ej }X ±iej } y)/y/2 are the circular components. 
The cross sections a a depend on u>, but not on j or 9 B - 

The total scattering opacity from mode j is = + 
k%, and the total extinction opacity 

In the diffusio n approximation (whose accuracy was 
studied, e.g., bv lShibanov fc Zavlinl 119951) . the effective 
opacity is 

nf = (cos 2 0/kJ + sin 2 d/nf)' 1 , (29) 

where d is the angle between B and the intensity gradi- 
ent, 

1 3 r cos 2 ^ 13 rwaPBg 

1 = 1 I a \ sm0 Bd9 B , — = -/ , r d9 B - 

(30) 



The effective opacity for the nonpolarized radiation is 
K ° s = 2/[(nf)- 1 + { K f)- 1 }. 

In a partially ionized atmosphere, the opacity is con- 
tributed by electrons, ions, and bound species. The scat- 
tering cross section includes contributions from the elec- 
trons and protons: <7® = a^ c + a^ p (the Rayleigh scat- 
tering by atoms can be important only at lower photon 
energies than considered in this paper). The absorption 
cross section c„ includes contributions from absorption 
by plasma electrons and protons (free- free transitions due 
to the electron-proton collisions, cr^, and proton-proton 
collisions, c^ p ), transitions between discrete states of an 
atom (bound-bound absorption, <7„ b ) and photoioniza- 
tion (bound-free, ffjj ). So, for the hydrogen atmosphere, 
we can write 

< = *h(C + a h J) + (1 - x H ) {al + <P), (31) 

where is the number fraction of atoms, which will be 
evaluated in the following section. 

3. EQUATION OF STATE 

3.1. Calculation of Tables 

Our treatment of the ionization equilibrium and EOS 
of hydrogen in strong magnetic fields is based on the the- 
ory developed in Paper I and briefly exposed in Sect. 12.21 
Since our free-energy model is computationally expen- 
sive, it is not possible to use the EOS code "on line" in 
any practical application. The alternative is to tabulate 
thermodynamic quantities covering the density, tempera- 
ture, and magnetic field domain of interest and to rely on 
an interpolation procedure in the table. Here we present 
EOS tables which cover a range of p, T, and B appro- 
priate for most typical neutron stars, such as isolated 
pulsars. 

As discussed in Paper I, our model becomes less re- 
liable at relatively low T and high p, particularly be- 
cause of formation of molecules and chains H„, which are 
treated in an approximate manner. In this domain, the 
partial number fractions and thermodynamic quantities 
are strongly model-dependent. However, this domain of 
uncertainty is unimportant for modeling of not too cold 
neutron stars, because the temperature grows inside the 
stellar envelope. For this reason, we have chosen as an 
input parameter the "astrophysical density parameter" 
R = Po/Tq, whic h is customary in the stellar opacity 
calcu lations (e.g., Ilglesias fc Rogers! Il996t ISeaton et al.l 
I1994I ). and restricted the calculation to R < 4 x 10 3 . 

3.2. Ionization Equilibrium 

Our detailed thermodynamic model shows that a 
strong magnetic field generally increases the fraction of 
bound species. In Fig. [I] ionization equilibrium curves 
at B = 10 12 G and 10 13 ' 5 G are compared with the 
case of B = 0. The latter case is treated in the frame- 
work of the nonmagnetic free ener gy model (Sect. II of 
Paper I), w hich is a variant of the Saumon fc Chabrierl 
([19911 H9921 ) model. In all cases, the excited atoms con- 
tribute significantly at low p. In the strong magnetic 
field, the population of decentered atoms is also signif- 
icant at low p. At higher density, the excluded-volume 
effect eliminates the excited and decentered atoms. At 
these high densities, the plasma species strongly inter- 
act, which leads to appearance of a significant fraction 
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density at T = 5 X 10 5 K and B = (dashed lines), 10 12 G (solid 
lines), and 10 13 5 G (dot-dashed lines). 



of clusters. Such clusters contribute to the EOS sim- 
ilarly to the atoms, lowering the pressure, but their 
radiation-absorption properties are clearly different from 
those of an isolated atom. Therefore they should be 
excluded from ih in Eq. (|31[) . Analogously, at low p 
we should not include i n xh the high l y exci ted states 
that do not satisfy the llnglis fc Teilerl (| 19391) criterion 
of spectral line merging, being strongly perturbed by 
plasma microfields. Such states form the so called optical 
pseud o-continuum (e.g., iDappen. Anderson. &; Mihalad 
1987). This distinction between the "thermodynamic" 
and "optical" neutral fraction is inevitable in the chem- 
ical picture of a plasma at high densities (see, e.g., 
iPotekhinl Il996bl for a discussion). We discriminate the 
atoms which keep their identity from the "dissolved" 
states (i.e., strongly perturbed by the plasma environ- 
ment) using the occupation probability formalism. At 
every s, v, and K±, we calculate the "optical" occupa- 
tion probability w° s (K±), replacing the Inglis-Teller cri- 
terion by an approxima te criterion based on the a verage 
atomic size [Eq. (14) of iPavlov fc Potekhinl (|1995l )]. The 
fraction of weakly perturbed atoms, which contribute to 
the bound-bound and bound-free opacities, constitutes 
a fraction w° s (K±)/w t l/s (K±) < 1 of the total number 
of atoms. Here, w^ s (K±) is the "thermodynamic" occu- 
pation probability derived from the free energy (Paper 
I). Heavy lines in Fig. [T] show the neutral fraction of 
the weakly perturbed atoms in their ground state, which 
contribute to the opacities as isolated atoms, whereas 
the light lines show the total fraction of protons bound 
in atoms or clusters. 

According to our model, at relatively low T, pres- 
sure ionization proceeds via a first-order phase transi- 
tion. This "plasma phase transition" occurs at temper- 
ature below T c » 3 x 10 5 B^ 9 K at densities around 
p c w 143i?i2 18 g cm' 3 (Paper I). In general, the validity 
of the free-energy models in the framework of the chem- 
ical picture of plasmas is questionable near the plasma 
phase transition domain. However, the T c and p c val- 



ues correspond to log 10 R c w 3.7, which is beyond the 
upper R limit for our tables, so that the plasma phase 
transition is not crossed along the tabulated isotherms. 

Figure [2] shows the domains of partial ionization in the 
T—p plane at three values of B. With increasing B, the 
domains where the atomic fraction is above a specified 
level expand significantly. For instance, at B — 5x 10 12 G 
the domain where xh > 0.01 extends to T — 10 6 K. Such 
amount of atoms can give an important contribution to 
radiative opacities. 

Our tables provide values of xh , as well as the fractions 
of ground-state atoms, molecules, and clusters at every 
R, T, and B entry. 



3.3. Thermodynamic Functions 

Figure [3] shows pressure P along two isotherms for the 
same field strengths as in Fig. [TJ The pressure varies over 
many orders of magnitude in the shown density range. 
Therefore, in order to make the discussed effects more 
visible, we plot in Fig. [3] the ratio of P to nok B T, the 
pressure of an ideal monatomic hydrogen gas at the same 
p and T. 

At different field strengths (including B = 0), the de- 
viations from the ideal gas behavior are qualitatively 
the same, but quantitatively different. At very low p, 
we have nearly fully ionized, almost ideal nondegener- 
ate gas of electrons and protons, so that P/n^k^T w 2. 
With increasing density, atomic recombination proceeds 
according to the Saha equation [in the strong magnetic 
field, the modified Saha equation is given by Eq. (54) of 
Paper I]. Therefore, the ratio P/nofceT decreases. At 
high density, however, the atoms become pressure ion- 
ized: in this region the increase of P due to the in- 
creased number of free electrons and protons competes 
with a negative nonideal contribution, which is mainly 
due to the Coulomb term in the free energy. The 
dot-dashed curve bends down near the high-p edge of 
the left panel of Fig. [3j because of enhancement of the 
nonideal contribution; it is a precursor of the plasma 
phase transition, where our model becomes inapplica- 
ble. Finally, at still higher densities, the electrons be- 
come degenerate and raise the pressure far above n^k^T. 
The dotted lines in Fig. [3J show the pressure of an ideal 
electron-proton gas; their upward bending marks the on- 
set of electron degeneracy. At B = 0, it occurs at 

p > m H {2m c k B T) 3 / 2 /(3TT 2 h 3 ) » 6T 6 3/2 gem" 3 , but in 
the strongly quantizing field the electrons become degen- 
erate at p > (mn/TT 2 ^) (m Q k B T '/2) 1 / 2 « 613 B 12 V% 
g cm -3 . 

Since the strong magnetic field enhances atomic re- 
combination and delays pressure ionization and electron 
degeneracy at high densities, the discussed features of 
the P/n$k B T curves become more prominent and shift 
to higher p with increasing B, as seen in Fig. [3l 

Along with the pressure, our EOS tables contain in- 
ternal energy U, entropy S, specific heat Cy, and the 
logarithmic derivatives of pressure xt — (d\n P/dlnT)v 
and Xp = ~ (d\n P/d\nV)T- Other second-order quanti- 
ties can be calculated us ing the Maxwell relations (e.g., 
lLandau fc Lifshit9 '.1993). For example, the heat capac- 
ity at constant pressure, Cp, and adiabatic gradient, 
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(<91ogT/<91ogP)s, are given by relations 



Cp = C\ 



PV 



' ad 



XT 



X 2 T + XpC v T/(PV) 



(32) 

Figure |4] shows V ac i at different values of B. At low 
density, the magnetic field increases the adiabatic gradi- 
ent, thus stabilizing the matter against convection. This 
thermodynamic effect is addition al to the hydrom a gnetic 
stab ilization considered, e.g., by iCha ndrase kharl ()196lD 
and iMiralles et abl (|1997| ). However, at higher densities, 
p ~ (1 — 100) g cm' 3 , there is a significant decrease of 
V a d due to the partial recombination of H atoms. The 
adiabatic gradient increases again at still higher densi- 
ties, where the plasma is fully pressure-ionized. 

4. CROSS SECTIONS 

4.1. Scattering 

The scattering cross sections by the electrons under 
the conditions typical for photospheres of the neutron 
stars with strong magnetic fields were thoroughly studied 



in the past (e .g., iKaminker et al.l fl982: Meszaroslll992t 
IVentural I1979L and references therein). The cross sec- 
tion cr a j_\ exhibits a resonance at the electron cyclotron 
frequency uj cc . Outside of a narrow (about the Doppler 
width) frequency range around uj ce , the cross sections are 



■ 0t, 



(33) 



where ctt = (87r/3)(e 2 /m e c 2 ) 2 is the nonmagnetic Thom- 
son cross section, and u e>a is an effective damping fre- 
quency given by Eq. (|38|) below. 

The scattering by the protons, which is completely neg- 
ligible in a nonmagnetized plasma, becomes important in 
the strong magnetic field, because aj^ exhibits a reso- 
nance at Lu cp . It is the same as the electron cyclotron 
resonance but a different mass and opposite charge of 
the particle: 



T s,p 



or- 



(34) 
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Fig. 4. — Adiabatic temperature gradient V a d at T = 5 X 10 5 
K, for different field strengths (shown by different line styles): 
log 10 S/G = 12.0, 12.3, 12.6, 12.9, 13.2, and 13.5. Triangles show 
the zero-field case. 



The damping frequency v v a will be derived below [Eq. 

We neglect the Doppler broadening of these reso- 
nances. Within the thermal width of the cyclotron reso- 
nanc e, the treatment of ra diation scattering is nontrivial 
(e.g. JVentura et aT]|1985l and references therein). How- 
ever, at T < 10 7 K, the Doppler width - uj ^JT/T, is 
smaller than the frequency resolution of our opacity ta- 
bles (chosen to be Alog 10 uj = 0.02). 

4.2. Absorption by Atoms 

Oscillator strengths for the H atom which rests 
in a strong magne tic field were calculated, e.g., by 
iForster et al.l (|1984f ). Bound-bound transitions of the H 
atom movi ng arbitrarily in a s t rong magnetic field were 
studied bv iPavlov fe Potekhinl (jl995). The modification 
of the binding energies due to the atomic motion (Sect. 
12. ip leads to a dramatic "magnetic broadening" of the 
spectral lines averaged over all states of motion, which 
exceeds by orders of magnitude the usual Doppler broad- 
ening. Thus the spectral profile of the bound-bound 
opacities becomes continuous in a wide frequency range, 
resembling a reversed bound-free profile. Our calculation 
of the bound-bound abs orption cross sectio n s reli es on 
the theory presented by IPavlov fc Potekhinl (|1995f) and 
employs fitting formulae for the binding energies, oscil- 
lator strengths, and electron collision widths derived by 
iPotekhinl (p98h . 

Photoionization cross sections of the nonmoving H 
atom in a magnetic field were calculated by many au- 
thors (e.g., lPotekhin. Pavlov, fc Vent ura 1997, and refer- 
ences therein) . Photoionization of the H atom in a strong 
magnetic field with allowance for motion was studied, 
using different mod ifications of the adiabatic appro xima- 
tion (Sect. I2TT1) bv iBezchastnov fc Potekhinl (| 19941) and 
Kopidakis, Ventura, & Hcrold (1996). A complete nu- 
merical treatment beyo nd the adiabatic approxim ation 
has been developed by iPotekhin fc Pavlovl (119971), who 
adapted the i?-matrix formalism ( Wiener fc Eisenbudl 



1947) to the case under study. They showed that none 
of the versions of the adiabatic approximation can pro- 
vide accurate crjj for all values of a and K±, partic- 
ularly because the continuum-channel coupling strongly 
affects o± f at sufficiently large K± . Here we use the com- 
plete numerical treatment. Since it is computationally 
involved, we use an interpolation across a precalculated 
set of tables. For each of the three basic polarizations, 
we have calculated a^ t i (uj,K±,B) on a predefined grid, 
with log 10 huj [eV] ranging from 1.0 to 4.5 with step 0.02, 
\og w K± [a.u.] ranging from 1 to 3 with step 0.1, and 
log 10 B[G] ranging from 11.9 to 13.5 with step 0.1. At 
the low-K± end of the grid, K± = 10 a.u., the atomic 
properties are virtually the same as at K = 0; beyond 
the upper limit, K± > 10 3 a.u., the contribution of the 
bound-free transitions to the total opacities is negligi- 
ble. The grid is sufficiently fine for calculation of atmo- 
sphere models. However, the step 0.02 in log 10 w does 
not allow us to resolve the narrow Beutler-Fano-type 
resonances which appear due to autoionizing states in 
the vicin ity of photoionization thresholds of partial cros s 
sections (jPotekhin fc Pavlovll997t IPotekhin etal]|1997l ). 
Whenever such a resonance occurs near a grid point, 
it produces a spurious outlier on the otherwise smooth 
cr„ f (uj, K±) dependence. We filter out such outliers by 
smoothing cr^ f as a function of uj at every K±, using 
the 3-point median filter. Since the grid does not allow 
us to resolve the photoionization threshold accurately, 
the threshold frequency wth is determined independently 
for every K \ , u sing analytic fits to the binding energies 
(|Potekhinlll998[) . 

In addition to the bound-bound and bound-free atomic 
transitions, in a plasma environment there are transi- 
tions from bound states to the highly perturbed atomic 
states discussed in Sect. 13.21 These perturbed levels 
effectively dissolve and merge in a pseudo-continuum, 
which lies below the photoionization threshold. In or- 
der to take into account the radiative transitions into 
this pseudo-continuum, we employ a below-threshold 
extrapolation as de s cribed for the zero-f i eld ca se by 
[ Dappen etlrtlti 19871) . IStehle fc Jacquemotj (|1993f ). and 
ISeaton et alj (|1994| ). Below w t h, the effective "bound- 
quasi-free" photoabsorption cross section due to the dis- 
solved lines is 



(uj < UJ th ) 



27re 2 < 



(35) 



m c c W( " duj 
where w° and are the optical occupation probabili- 
ties of the initial and final states, respectively, /i/ jQ is 
the corresponding oscillator strength, and d^//dw is the 
number of fi nal states per unit frequen cy interval. By 
analogy with lStehle fc Jacquemotl (|1993f ). we interpolate 
w°f as function of frequency and set w° f = w° at w smaller 



than the lowest allowed transition frequency. Taking into 
account that o~\ (uj < u>th) is a smooth continuation of 
a^(u> > w t h) (|Da,ppen et al.lll987l ). we write 



(36) 

where <7„ > extr (A"^, cu) is a power-law extrapolation of 
a^(K±,uj) at uj < u t h- 

Unlike the B = case, in our case fif, a an d de- 
pend on polarization. For photoabsorption by an atom 



Atmospheres of magnetic neutron stars 



963 



in the ground state, /i/,±i 7^ only for even upper states, 
whereas /i/,o 7^ only for odd states. We take into con- 
sideration only the appropriate states while calculating 



w° f in Eq. (36J). 

The if^-dependent cross sections are averaged over the 
distribution of atoms over with statistical weights 
w°(K ± ) ex p\-e l (K±)/k B T] 1 as in | Pavioy_ &: Potekhinl 
(fl99fih and lPotekhin fc Pavlovi (fl99l . 

4.3. Free-Free Absorption 

In the class ical cold plasma approximation (e.g., 
lGinzburelll970l ), the free-free absorption by electrons is 



1 



Aire 2 



(uj + auj cc ) 2 + i/f m c c 



(37) 



where 1/? is an effective frequency of electron-proton col- 
lisions which lead to absorption of photons. Broadening 
of the electron cyclotron resonance in Eqs. (f3"3"| and (|3T|) 
is determined by the sum of the effective frequencies for 
absorption and scattering, 

(38) 



vl + vt 



where 



2 e 2 



3 m e c 3 



(39) 



is the natural (radiative) width of the resonance. The 
damping frequency given by Eq. (|38p ensures the correct 
value of the cyclotron abso rption cross se ction integrated 
across the resonance fe.g.. lVenturaHl979h : 



"cc+A" 47r 2 e 2 
[er^(w) + a s 1 (uj)] duo — , 

' cc — Aui m c C 



(40) 



The values of z/f are provided b y quantum-mechanical 



with i/ 6ja <C Aw <C oj c 

3 provided byj 

calculations. It is customary fe.g.. lArmstrong fc NichoTlsl 
|1972|) to express cr ff through the thermally averaged 
Gaunt factors g, or equivalently, Coulomb logarithms 
A = (Tr/y/E)g. Taking into account Eq. (|37|) . we can 
write 

- (l-e-»)A« 



3 V m c k^T huj 



Hco 



(41) 



The factor (1 — e ") allows for the induced radiation. 

4.3.1. Infinite Proton Mass Approximation 

In the zero-field case, the electron free-free absorption 
rate can be calculated assuming the electron scattering 
off a fixed Coulomb potential. In this case, allowance 
for the finite ion mass consists in replacing m e by the 
reduced mass fj,. The Born approximation yields the 
well-known formula for the cross section of free-free pho- 
toabsorption by an electron having an initial momen- 
tum pj and final mome ntum pj — [p 2 + 2m c huo) 1 ^ 2 (e.g., 
iBethe fc Salpeterlll957ft . 



<r s (Pi,uj) = 



I6n 2 



3 m e choj 3 p, 



In 



Pf + Pi 



Pf-Pi 



(42) 



whose averaging over the Maxwell distribution gives the 
classical Coulomb logarithm 



a: 



' 1,2 K {u/2), 



(43) 



Kq(u/2) being the modified Bessel function. iHummerl 
(|1988h calculated g e using accurate non-Born quantum- 
mechanical results bv lKarzas fc Latter! (|1961[ ) and fitted 
it by a Pade formula. The nonmagnetic Gaunt factor 
is applicable if the magnetic field is nonquantizing - 
i.e., if /3 C < 1, where /3 e is given by Eq. (fTJ). In the 
quantizing magnetic fields, the Coulomb logarithm was 
evaluated in the Born approximation by several au thors 
(|Meszaroslll99l iNagell fl980l : iPavlov fc Panovjll976l ). In 
this approximation, A*, which is generally a function 
of B, T, and u>, depends only on the two dimensionless 
arguments, u and (3c 1 : 



A ff = - e 

4 e 



l/2 ^ 



n— — 00 



l%((3c,u,y) dy, 



where 



At 



A 



:(/3 °' W,y) C [(y + + C)sinh(/3 c /2)]l« 
o _ x n Ki{x n ) ±1 _ y + 6 + \n\( 



(44a) 



(44b) 



V + &/4 



C= \/l + 2% + y 2 , 



C 2 

■ exp(-/3 c ) 



1 - exp(-/3 c ) 



K (x n ){Uc) 
(44d) 



x n = \u- n/3 e \ y/0.25 + y//3 . (44c) 

Equation (|4"4"1) has been derived assuming that an elec- 
tron scatters off a fixed Coulomb center. Actually the 
protons are moving and can absorb radiation during col- 
lisions. Although this process is negligible at B = Q, it 
may be important at uj near or below uj cp . The previous 
authors (THo fc Lail[200lt [C^[20f)ll : IPavlov et al.lll995l : 
IZane et aLll2000L |2001[ ) supplemented erf by a cross sec- 
tion of "ion free- free absorption," crf' p . For hydrogen, 
taking into account Eqs. (|33[) and their formulae 

can be written as 

< T a = (J a + a a ) a / a = u otl a a • ( 45 ) 

We find, however, that Eq. (|45|) is erroneous. 



4.3.2. Absorption in Proton Collisions 

There are two effects of the finite proton mass on the 
absorption: first, the absorption can occur in proton- 
proton collisions, and second, the absorption in the 
electron-proton collisions is modified because of the pro- 
ton motion. Let us start with the first process. By anal- 
ogy with Eq. (|37p . we write the cross section as 



1 



(uj - auj cp ) 2 + v\ 



p.a 



47re 2 



m n c 



(46) 



where v pp is the effective frequency to b e determined. 

In the classical picture (|Ginzburglll970l ). an ion- ion col- 
lision does not change the total electric current and hence 
does not cause dissipation. Therefore, it does not con- 
tribute to the damping of radiation. In quantum mechan- 
ics, this corresponds to vanishing dipole matrix element 
for the absorption. In Appendix fA~l we evaluate the non- 
vanishing quadrupole term under the condition that the 

1 The set of equations i144l) is equivalent to Eq. (27) of 
IPavlov fc Panovl [IT976T ), but in Eq. iHIcll for A^ 1 we have restored 
power 2 of lacking in Pavlov & Panov (1976) apparently due to 
a misprint. 
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Fig. 5. — Coulomb logarithms for photoabsorption in collisions of 
nonrelativistic protons (A pp , solid line), collisions of identical but 
distinguishable particles (dashed line), and electron-proton scat- 
tering (AS, dotted line). 



magnetic field does not quantize proton motion, that is 
/3 P < 1, and obtain 



,pp 



256 
3 



m-pk^T fku m p c 2 



(l-e' u )A 



pp. 



(47) 



where A pp is an appropriate Coulomb logarithm. In Fig. 
O App is plotted by the solid line. For comparison, we 
also show A for distinguishable particles [corresponding 
to Eq. (|A4[k dashed line], and the classical Coulomb log- 
arithm Afj [Eq. ([HD; dotted line]. The fit 



App « 0.6 ln(22u" 



9 m" - 3 ) +0.4^ (48) 

accurately reproduces A pp at small and large u and has 
a maximum error within 1.5% at intermediate u. 

We see that a pp differs from er* ,p [defined by Eq. 
(145]) with use of Eqs. ([37]) and Ipflj) ] by a factor of 
32 v/2m p /m c (k B T/m p c 2 ) A pp /A% ~ T/10 9 K. 

One should remember that Eqs. (|46 |) - (|48l) are obtained 
in the nonrelativistic approximation. Therefore they do 
not take into account spin-flip processes and are inap- 
plicable at very high T or high uj, where the relativistic 
corrections can be important. 

4.3.3. Electron Free-Free Process with Allowance for 
Finite Proton Mass 

Let us write the cross section of photoabsorption due 
to the electron-proton collisions in the form of combined 
Eqs. (|37|) and (|4ip . neglecting broadening: 



2 9 / 2 tt 3 / 2 n e e 6 



1 



A? 



3me /2 (fc B T) 1/2 c?iw (uj + auj cc ) 2 a 



(49) 



The superscript 'e' indicates that the electron (not pro- 
ton) cyclotron resonance has been separated off A a . 
Since the colliding electron and proton are treated on 
equal footing, we anticipate that thus defined normalized 
cross section A^ will reveal a resonant peak at w ~ cj cp . 

The initial and final states of the interacting electron 
and proton are just continuum states of the H atom. An 
accurate treatment of these states would imply a solution 
of the coupled-channel equations (fl5|) and calculation of 



the i?-matrix, as we did for the bound-free process. How- 
ever, we will restrict to the first Born approximation. In 
this approximation, A° is given by Eqs. (|B20[) - (|B22[) de- 
rived in Appendix [Bj For the longitudinal polarization 
(a = 0), a calculation by these equations well reproduces 
the Coulomb logarithm (|44|) obtained in the infinite pro- 
ton mass approximation. However, for the two circular 
polarizations the result is different; it is shown in Fig. [6] 
by the solid lines. As expected, we see a remarkable pro- 
ton cyclotron resonance at a = +1, which is due to the 
denominator uj — uj cp in the last terms of Eq. (|B22c[) . In 
addition, for both circular polarizations there are smaller 
spikes at higher proton-cyclotron harmonics, arising from 
the l ogar ithmic singularities oi v nsn i 8 {p, k) at x — ► [see 
Eq. (|B8[) ]. Apart from these spikes, A^ is accurately de- 
scribed by the formula A±! w A± x ui 2 / (u — aLu cp ) 2 , where 
A^ is given by Eq. (|4"4"|) (dot-dashed lines in the figure). 
We note that the factor u) 2 j(to — ao; cp ) 2 naturally ap- 
pears in the classical plasma model with allowance for 
the ion motion. The classical model also helps to re- 
store the damping factors neglected in Eq. (|49[) . Let v e 



and v p be the electron and proton damping frequencies 
due to processes other than the electron-ion collisions. 
For the processes considered above (Thomson scatter- 
ing and proton-proton collisions 2 ), we have v a = and 
v p = v p + i/ pp , where v p — 2e 2 uj 2 / (3 m p c 3 ) is the nat- 
ural width of the proton cyclotron resonance. Averaged 
Newtonian equations of motion for the electrons and pro- 
tons in the magnetic and rad iation fields give the com- 
plex permittivity tensor e (cf. lGinzburglll970l §10). Ne- 
glecting ^/m e /m p compared to unity and assuming that 
v p <C v c <C uj, we have 



uj 



lm(e xx + iae xy ) = — 

UJ 



^ cp (^cc + Cfw) v p + uj (lu - auj cp ) v c 

I [ (lu + o.lu cc ) (lu - auj cc ) 

~ v a { v v + v c m c /m p ) -v p v c ] 2 
+ [lu v^, + (w + aui cc ) v p + (uj - au) cp ) v e ] 2 \, (50) 

where a — ±1, and Im means the imaginary part. Ne- 
glecting the tiny shift of the position of the cyclotron 
resonances caused by the damping, we can now write 



lu 2 47re 2 z/f 
(uj + auj cc ) 2 (uj ~ auj cp ) 2 + lu 2 i> 2 m c c 



(51) 



where 



v a = + (1 + auj cc /uj) Up + (1 - auj cp /uj) v c , (52) 

and is given by Eq. (|4Tj) . 

The accurate calculations according to Appendix 151 are 
reproduced by Eqs. (fSTj) and (|4"Tj) . if to multiply A^ in 
the latter equation by a correction factor of the order 
of unity. This quantum correction factor proves to be 

2 We do not consider the electron-electron collisions. At B = 0, 
they are known to be unimporta nt for the bremsstrahlu ng, ex- 
cept for relativistic energies (e.g., Bcthc & Salpctcr 1957). The 
magnetic field does not change this conclusion, because the res- 
onance, that appears at lu = LU ce for this process, merges in the 
more powerful classical cyclotron resonance for the usual free-free 
and Thomson processes. 
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Fig. 6. — Comparison of acc urat e and approximate normalized 
free-free cross sections [Eq. | |49| |1 for two circular polarizations 
a = ±1 at 3 e = 10 4 . Solid line: Eqs. (fB19t fB22l ; dot-dashed 
line: Eq. ( 1511 : short-dashed line: customary approximation, Eq. 
(|45T : long-dashed line: A*, Eq. Bit : dotted line: A^, Eq. | [43t . 
Arrows indicate the proton cyclotron harmonics. 



the same for a — +1 and a = — 1. Thus, the two effec- 
tive collision frequencies (longitudinal and transverse 



+i 



s^i) provide the three 



Near the electron cyclotron resonance (to ~ uj cq , a = 
— 1), the damping frequency v a approximately repro- 
duces v e ^ a in Eq. (|38p . which ensures the condition 
(|40|) . Near the proton cyclotron resonance (ui w ui cp , 
a = +1), the effective damping frequency is v a ui/uj cc m 



Y m e /m p 



where 



yff.P 



(53) 



Equation (|53p is consistent with the requirement of os- 
cillator strength conservation for the proton cyclotron 
resonance in polarization a = +1, fully analogous to Eq. 




3.5 4 
log hw (eV) 



4.5 



Fig. 7. — Monochromatic radiative opacities for the ordinary 
(dashed lines) and extraordinary (solid lines) polarization modes 
in a plasma at p = 500 g cm -3 , T = 5 X 10 6 K, and B = 5x 10 14 
G, for 9g = 0, 1 0°, and 90°. Dotted lines show the opacities 
according to Eq. 11451 . Here the opacity is calculated assuming 
complete ionization. 



From the relations 



= 1.78 x 1CT 4 T 6 



A 



,ff, P 



pp 
A ff ' 



3.6po 1 — e 



T, 



5/2 



A 



ppj 



(54) 

we see that the proton-proton collisions can be safely 
neglected at any T and p typical of outer envelopes of 
the neutron stars. 

For comparison, Fig. [5] also shows the nonmagnetic 
Coulomb logarithm A^ (dotted line), A^ given by Eq. 
(|44|) (long dashes), which neglects the finite proton mass, 
and A^ which would correspond to the calculation of er^ 
according to the traditional Eq. (I45|) (short dashes). It 
is easy to see the difference of Eqs. (l4T?|) and (|5ip from 
Eq. (|4"S"1) . The customary recipe (|4"5)) miss es the interfer- 
ence of the first two terms in each of Eqs. (|B22b[) . (|B22c[) 
(related to transitions with changing Landau number n) 
with the terms in square brackets (related to transitions 
which keep n constant). At w < w cp the latter terms 
tend to compensate the former ones. A suppression fac- 
tor ~ (lu/lu cp ) 2 , which follows from Eq. (|5"Tj) at uj <C w cp , 
is brought about by this interference. From the classical 
physics point of view, it may be explained as follows: very 
slow (ui <C 0J C p) oscillations of the radiation electric field, 
perpendicular to the constant magnetic field, make both 
particles, electron and proton, to drift adiabatically with 
the velocity (c/B 2 )E x B, so that in the drifting frame 
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of reference they do not "feel" the electric field of the 
electromagnetic wave. The suppression of the free-free 
cross sections takes place for both circular polarizations. 

As well known, the Born approximation is accurate 
only at Hlo much larger than the binding energies. In 
order to partly correct A^ beyond this approximation 
and to recover the non-Born Gaunt f actor q s at B = 0, 
we multiply A^ by the ratio of g s (Hu mmer] I1988D to 
SBom = (V3/7r)A* 

5. OPACITIES 

5.1. Fully Ionized Atmosphere 

For a fully ionized atmosphere model, the monochro- 
matic opacities calculated in the framework of the theory 
outlined in Sect. l2.3l do not differ m uch from the opacities 
used in the previous models (e.g.. IShibanov et"aLlll992f l 
at ijj > u! cp . However, the improved treatment of the free- 
free contribution results in a considerable modification of 
the opacities at u < lu cp . Figure [7] shows monochromatic 
opacities for two polarization modes, j = 1 (solid lines) 
and j = 2 (dashed lines) for radiation propagating at 
three angles 9 b with respect to the field lines. The atmo- 
sphere parameters chosen for this figure can be expected 
near the bottom of a magnetar photosphere: p = 500 
g cm" 3 , T = 5 x 10 6 K, and B = 5 x 10 14 G. An exten- 
sion of our EOS for the s uperstrong magnetic fields (work 
curre ntly in progress — IChabrier, Douchin. &: Potekhinl 
i2002T) indicates that the model of a fully ionized at- 
mosphere is adequate for this relatively high value of 
p. At these parameters, the proton cyclotron resonance 
at tuv = 3.15 keV is quite prominent. At small 9b, 
the opacities of the two normal modes cross each other 
at tu sa uy p , which i s a well known phenomenon (e.g., 
IShibanov et al.l 119921 ) . Another mode crossing, which 
occurs at Huj sa 7.6 keV, is due to the vacuum reso- 
nance (e.g., Pavlov fc Gnedinlfl98l IShibanov et al.f l992: 
iSoffel et "all 1 19831 ). Near the crossing points the modes 
may become completely nonorthogonal ("collapse") at 
certain angles, so that their designation is ambiguous 
(|Soffel et al.l[l98l . 

The dotted curves in Fig. [7] are obtained using Eq. 
([45)) . We see that this traditional calculation strongly 



overestimates K\ at large 9b, and overestimates both k\ 
and K2 at 9b = 0, if fka is small enough. 

Figure [8] shows the effective opacities n eS for the dif- 
fusion of nonpolarized radiation along (||) and across 
(_L) magnetic field (Sect. 12. 3| . for the same plasma pa- 
rameters as in Fig. [7] The peak at log hui/eV = 3.5 
is due to the ion cyclotron resonance, and the one at 
log huj/eV — 3.9 due to the vacuum resonance. The 
barely visible intermediate spike at log fku/eV = 3.8 is 
the quantum resonance of the Coulomb logarithm at the 
doubled proton-cyclotron frequency (cf. Fig. [S]). 

We see that the improvement of the free-free cross sec- 
tion discussed in Sect. 14.3151 is important for the effective 
opacities shown in Fig. [5] at u> < u> cp . At u) < 0.3 uj cp , 
the difference exceeds one order of magnitude. More- 
over, it has an impact on the Rosseland mean opacities, 
as discussed below. 

5.2. Partially Ionized Atmosphere 

As follows from Sect. 13.21 the amount of neutral hy- 
drogen in neutron-star photospheres can be significant 
at T < 10 6 . For example, at p = 0.1 g cm" 3 , T = 10 5 5 
K, and B = 2.35 x 10 12 G, 12% of protons are bound in 
the ground-state H atoms. Monochromatic opacities for 
this case are shown in Fig. [5] for three basic polarizations. 
The contribution of the fully ionized plasma component 
is shown by dot-dashed lines, while dotted and dashed 
lines show bound-bound and bound-free contributions, 
respectively. The total opacities are plotted by the solid 
lines. We see that the bound-free contribution is im- 
portant for any polarization, whereas the bound-bound 
opacity is important for a = ±1 but unimportant for 
the longitudinal polarization (a = 0). This is because 
the dipole selection rule forbids radiative transitions be- 
tween different tightly bound states with absorption of 
a photon polarized along B. Transitions to the odd hy- 
drogenlike states {y = 1,3, . . .) are allowed, but the cor- 
responding occupation probabilities are small, so that 
these upper levels are effectively merged into the contin- 
uum. The absorption peak at log 10 hui/eV ~ 1.84 on the 
left panel corresponds to the transition from the ground 
state to the state with s = 1 for the centered atoms. It 
would be a narrow spectral line without atomic motion. 
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Fig. 9. — Monochromatic opacities for the basic polarizations (a = +1, left panel; a = —1, middle panel; a = 0, right panel) in a typical 
partially ionized neutron-star atmosphere at p = 0.1 g cm" " 3 , T = 10 5 - 5 K, and B = 2.35 X 10 12 G. Dot-dashed lines: opacities of fully 
ionized component; dashed lines: bound-free opacities; dotted lines: bound-bound opacities; solid lines: total opacities. 
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Fig. 10. — Monochromatic opacities for the same basic polarizations as in Fig. M at T 
-3 n „ „™-3 and 5 q cm -3 ] and B _ 1Q 13 Q 
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lines), p = 0.05 g cm 3 , 2 g cm 



2 2.5 3 3.5 
log hu (eV) 

10 5 K (solid lines) and 2 X 10 6 K (dashed 



On the middle panel (a = —1), there would be no sig- 
nificant bound-bound absorption at all, were the motion 
neglected. However, the thermal motion of atoms across 
the field drastically modifies the spectrum. There ap- 
pears significant absorption for a = — 1. The bump at 
fojj m 100 eV for a = +1 is due to the transition to the 



second excited level (s = 2), which would be negligible 
for nonmoving atoms. The magnetic broadening men- 
tioned in Sect. 14.21 smears the photoionization edges at 
huj > 100 eV and extends the bound-bound absorption 
frequency range down to ~ 10 eV for any polarization. 
The spikes near this low to are explained by the K±- 



968 




O 
CO 

a 

o 

cu 
^> 
ID 
a 
cu 



em 
o 




-4 -3 -2 -1 
log p (g/cm 3 ) 



4 -3 -2 -1 
log p (g/cm 3 ) 



Fig. 12. — Effective Rosseland mean opacities for diffusion of nonpolarized radiation along (solid lines) or across (dashed lines) magnetic 
field, at T = 10 5 5 K (left panel) and 10 6 5 K, for B = 0, 10 12 G, and 2 X 10 13 G. Dotted lines show the opacities of a fully ionized plasma 
(at each value of B ^ 0, the lower dotted curve corresponds to the transverse diffusion and the upper to the longitudinal one). 



depen dence of the transition energy (|Pavlov &: PotekhTnl 
I1995D . The spike at w w w cp (log 10 fiw/eV w 1.17) 
appears because the transitions between the decentered 
states whose quantum numbers s differ by one (As = 1) 
correspond to the energies « huj cp , almost independent 

of K± [it follows from Eq. (|14|). since e\ v (K±) is small 
at large K±\. Another resonance occurs at a frequency 
slightly below cj cp , which corresponds to the minimum 
distance between the states with As = 1 (the avoided 



crossings, cf. lPotekhirilll994| ). because the energy differ- 
ence has zero derivative with respect to K there. These 
resonances are smoothed by the electron impact broad- 
ening. 

In general, we see that partial ionization must be taken 
into account in the opacity calculations at these plasma 
parameters. 

Figure [TO] shows the total extinction coefficients pn a 
for different p and T at B = 10 13 G. The curves are 
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truncated from the left at w = uj p \. At p — 50 g cm -3 , 
the opacities are smooth functions of u, which reflects the 
fact that virtually all excited atomic levels are merged 
into continuum at these p and B. At p < 2 g cm" 3 , 
however, the curves clearly reveal the features due to the 
bound-bound and bound-free transitions. 

The opacities for the three basic polarizations, com- 
bined with appropriate components of the polarization 
vectors, provide the opacities in the two normal modes 
(Sect. l2~3")) . For example, Fig. [TT]shows the normal-mode 
opacities for p = 0.1 g cm" 3 and B = 10 13 G, at three 
values of 9 b and two values of T. At the lower T = 5x 10 5 
K, the features arising from the bound-bound and bound- 
free transitions are clearly visible at any 9b- 

5.3. Rosseland Mean Opacities 

Along with the thermodynamic functions and number 
fractions of species, our tables contain Rosseland mean 

effective opacities for longitudinal (/%) or transverse 
( k r) propagation of nonpolarize d radiation. They are 
calcu lated in a standard way (e.g.. lArmstrong fc Nicholisl 
Il972[ ) from the effective monochromatic opacities for the 
diffusion approximation, k" and ft 4- , defined in Sect. [2~51 
The improvement of the free- free cross section substan- 
tially affects «4 and k 4 , . For example, in the case of fully 
ionized plasma shown in Fig. [51 we obtain re^ = 1.8x 10~ 4 
cm 2 g _1 and k 4 , = 1.6 x 10~ 4 c m 2 g~ 4 , in reasonable 
agreem ent with the analytic fit in Potekhin fc Yakovlevl 
(pOOll ) (2.1 x 10" 4 cm 2 g" 1 and 1.7 x 10" 4 cm 2 g" 1 , re- 
spectively) . Since the latter fit did not take into account 
the ion cyclotron resonance, we conclude that this reso- 
nance is unimportant for the Rosseland opacities in the 
given example. Meanwhile, the traditional treatment of 
this resonance [Eq. |45)) ] yields effective opacities shown 
by the dotted lines in Fig. [8j whose Rosseland means are 
K;| = 1.3 x 10~ 3 cm 2 g" 1 and k 4 , = 1.0 x 10~ 3 cm 2 
g _1 , — that is, about six times larger than the accurate 
values. 

Figure [12] illustrates the density-dependence of 
the Rosseland opacities at two values of T for 
B = 10 12 G and 2 x 10 13 G. For comparison 
the nonmagnetic OPAL opacities (jlglesias fc Rogers! 
119961 |http : //www-phys . llnl . gov /Research/OPAL/) 
are showm 3 Dotted lines represent tne effective 
Rosseland opacities in the model of a fully ion- 
ized electron- proton plasm a, according to the fit of 
iPotekhin fc Yakovlevl (|200il ). 

As well known, the strong magnetic field makes the at- 
mosphere more transparent at given p and T, because of 
the presence of large w cc in denominators of Eqs. (|33|) and 
(j4T))) . At sufficiently large p, there is a good agreement 
between the opacity tables and the analytic fully-ionized 
plasma model. However, at p < (1 — 10) g cm -3 , there 
are large differences, which reflect the contribution of 
bound-bound and/or bound-free transitions in the effec- 
tive opacities. Remarkably, this difference is significant 
at T = lO 5 5 K even in the nonmagnetic case (the up- 
per curve on the left panel). As noted in Paper I, the 

3 As mentioned in Sect. IID of Paper I, our model of partially 
ionized hydrogen plasma at B = accurately reproduces the OPAL 
opacities. 



contribution of bound species increases with increasing 
B. 

In the case of T = 10 6 5 K (right panel of Fig. [T2"j) . 
the model of fully ionized plasma is quite accurate at 
B < 10 12 G. However, this is not the case at the higher 
field strength B — 2 x 10 13 G, where the contribution 
from bound species is again appreciable. 

6. DESCRIPTION OF THE TABLES 

The input parameters for our tables are B, T, and the 
astrophysical density parameter R — po/T^. At present, 
the tables are calculated for 11.9 < log 10 B/G < 13.5 
with step Alog 10 B = 0.1, 5.3 < log 10 T/K < 7.0 with 
step Alog 10 T = 0.05, and -7.4 < log 10 R < 3.6 with 
step Alog 10 i? = 0.2. 

The tables for different values of T and B have identical 
structure. An example is shown in Tabic [T] The first 
line contains log 10 T/K and Iog 10 i3/G. Each row then 
provides: 

1. log 10 i?; 

2. log 10 P, where P is the pressure in bar= 10 6 dyn 
cm -2 ; 

3. the dimensionless pressure parameter 
PV/ (Nok-BT), where Nq is the total number 
of protons (free and bound) in volume V; 

4. the dimensionless internal-energy parameter 
U/(N Q k B T); 

5. the dimensionless entropy parameter S / (Nobs); 

6. the reduced heat capacity Cy / (N kB); 

7. the logarithmic pressure derivative xt = 
(dlnP/dlnTV; 

8. the logarithmic pressure derivative Xp = 
(91nP/91np) T ; 

9. the atomic fraction xh, that is the total number of 
H atoms with non-dissolved energy levels, divided 
byiVo; 

10. the ground-state atomic fraction; 

11. the molecular fraction (the number of H2 molecules 
with non-destroyed levels, divided by No); 

12. the fraction of protons comprised in clusters and in 
strongly perturbed atoms and molecules; 

13. lo 

Sio K "i where k" is the effective Rosseland mean 
opacity for transport of nonpolarized radiation 
along magnetic field lines in the diffusion approxi- 
mation, in cm 2 g ; 

14. log 10 K 4- , where k 1 is analogous to k", but for dif- 
fusion of radiation in the direction perpendicular 
to B. 

We have also written a Fortran program for 
the cubic-polynomial interpolation of the tabulated 
data in the 3-parameter space of B, T, and 
p. The tables and the program are available at 
http : //www . iof f e . rssi . ru/astro/NSG/Hmagnet/ 
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Table 1. Sample EOS and opacity table. 







thermodynamic functions 




number fractions 




log(opacities) 


(1) 


(2) 


(3) (4) (5) (6) (7) (8) 


(9) 


(10) (11) 


(12) 


(13) (14) 



6.000 13.000 

-7.40 0.8173 1.999 1.256 59.737 2.378 1.001 1.000 3.77E-04 7.15E-05 0.00E+00 7.81E-03 
-7.20 1.0173 1.999 1.256 58.818 2.378 1.001 1.000 4.15E-04 8.30E-05 0.00E+00 8.85E-03 



-0.20 
0.00 
0.20 
0.40 
0.60 
0.80 



7.9984 
8.1949 
8.3905 
8.5853 
8.7793 
8.9721 



1.914 
1.898 
1.879 
1.857 
1.831 
1.801 



1.028 26.662 2.473 1.062 0.984 1.19E-02 8.10E-03 6.31E-11 9.91E-02 



0.977 
0.915 
0.841 
0.755 
0.654 



25.736 2.501 
24.807 2.536 
23.873 2.575 
22.935 2.616 
22.000 2.658 



1.074 0.980 

1.088 0.976 

1.105 0.971 

1.123 0.966 

1.142 0.963 



1.33E- 
1.45E- 
1.51E- 
1.45E- 
1.23E- 



9.09E- 
1.01E- 
1.08E- 
1.09E- 
9.97E- 



1.77E- 
4.26E- 
8.15E- 
1.09E- 
8.36E- 



10 
10 
10 
09 
10 



1.06E-01 
1.13E-01 
1.21E-01 
1.32E-01 
1.44E-01 



3.00 10.9665 1.122 

3.20 11.1018 0.967 

3.40 11.2407 0.840 

3.60 11.4056 0.775 



-1.917 11.401 2.765 

-2.457 10.364 2.836 

-3.068 9.334 2.900 

-3.736 8.301 2.946 



1.935 
2.276 
2.651 
2.906 



0.692 
0.670 
0.734 
0.975 



0.00E+00 0.00EH 
0.00E+00 O.OOE-j 
0.00E+00 O.OOE-j 
0.00E+00 O.OOE-i 



00 0.00E+00 1.50E-02 

00 0.00E+00 3.76E-03 

00 0.00E+00 4.11E-04 

00 0.00E+00 1.40E-05 



-4.896 
-4.872 

-1.242 
-1.091 
-0.954 
-0.835 
-0.738 
-0.662 

2.207 
3.083 
4.248 
5.782 



-4.896 
-4.872 

-1.447 
-1.295 
-1.158 
-1.039 
-0.941 
-0.866 

2.039 
2.904 
4.055 
5.574 



7. CONCLUSIONS 

We have calculated the EOS and radiative opacities 
of fully and partially ionized hydrogen plasmas in a wide 
range of densities, temperatures, and magnetic fields typ- 
ical for photospheres of the strongly magnetized neu- 
tron stars. The first- and second-order thermodynamic 
functions, non-ionized fractions, and effective Rosseland 
mean opacities are published in the electronic form. 

The opacities are calculated more accurately than in 
the previous publications. In particular, we take into ac- 
count suppression of the free-free absorption below the 
proton cyclotron frequency, which was overlooked pre- 
viously. This effect reduces the opacities of the ion- 
ized component of the plasma by orders of magnitude 
at photon energies Hco < 0.3 Huj cp ~ 0.02 Byi keV, which 
necessitates a revision of the prev iously published mod- 
els of X-ray spectra of magnetars (|Ho fc Laill200lL 120031 : 
IOzelll200lL 120031 : IZane et al.ll2001[ ). On the other hand, 
the bound-bound and bound-free absorption, neglected 
in the previous models of neutron-star atmospheres, in- 
crease the opacities by more than one order of magnitude 



at hoj ~ (0.1-3) keV in the outer atmosphere layers of 
the ordinary neutron stars with B ~ 10 12 -10 135 G and 
T < (1-3) x 10 6 K, which can also significantly affect the 
spectra. 

One can expect that the effect of the bound species 
on the EOS and opacities is as important for magnetars 
(despite their supposedly higher temperatures) as for the 
ordinary neutron stars. To check this, we need to ex- 
tend our model to hi gher B; preliminary high-B results 
(jChabrier et al.ll2002f) support this anticipation. 
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APPENDIX 

A. PHOTOABSORPTION DUE TO PROTON-PROTON COLLISIONS 

Th e general formula for the dif ferential cross section of absorption of radiation by a quantum-mechanical system is 
(e.g.. lArmstrong fc Nich"oTIs1ll972t ) 



da = ^- |e • (/|j cff |z)| 2 6(E f - Ei - hw)dv f , 



(Al) 



where i and / are the initial and final states of the system, dvf is the density of final states, e is the polarization 
vector, j e ff = e lk ~i' v j, j is the electric current operator, and fc 7 is the photon wave number. For two charged particles 
in a magnetic field, 

jeff = (e ik - ri qi wi/mi + e ik ^q 2 ir 2 /m 2 ), (A2) 

where qi and m t are the particle charge and mass (i = 1,2), and T?i is given by Eq. ([3]). Introducing, in the standard 
way, the center-of-mass (R, P) and relative (r, p) coordinates and momenta of two protons and using the gauge 
A(n) = (1/2) B x n, we get 



Jeff : 



Ak-,R e 



-B x R + Uk^ 



■r)(p 



4c 



B x r 



(A3) 



The first two terms do not contribute to the free-free absorption, because they do not contain the relative variables and, 
therefore, are decoupled from the Coulomb interaction. The remaining terms are similar to j c g relevant to absorption 
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of radiation by a particle with charge e/2 and mass m p , except for the factor fc 7 • r, which is small at hui <C m p c 2 . 
The absorption cross section can be written in the form of Eq. (|46|) . The nonquantizing magnetic field, j3 p <C 1, 
does not affect the effective collision frequency i/P p , which in this case does not depend on a. We have evaluated 
the proton free-free cros s section in the nonrelativistic Born approximation, using the technique of Fourier transforms 
(Bct he fc Salpeterlll957l . §77). For two distinguishable particles of equal mass m p and charge e, the cross section is 



256tt 2 



[ f PiPf , 3 P 2 +P 2 f 

rripchui 3 pi I (m p c) 2 10 (m p c) 2 



In 



Pi + Pi 



Pf - Pi 



(A4) 



Taking into account that the colliding protons are identical and have the spin 1/2, one should calculate the matrix 
element in Eq. (|A1[) for symmetric and antisymmetric final stat es and sum up the cross sections with the statistical 
weights 1/4 and 3/4, respectively (e.g.. lLandau fc Lifshitdfl976l ). This leads to the equation 



where 



1287T 2 



n D e u 



1/6 PiPf 



4 i 4 
Pi+Pf 



3 m p chuj 3 pi \ 5 (m p c) 2 (pf + p 2 f) 2 



c 



2 ^4 



pf+pj + 0.8 pjpj 



(pf+pjf 



hi 



Pf + Pt 



Pf - Pi 



(A5) 



(A6) 



In Eqs. (|A4[) and (|A6|) . p 2 = p 2 +m p huj, since the reduced mass equals m p /2. The Maxwell distribution for the relative 

momenta is .F pp (pi) = (4/y / 7r) (m p kBT)~ 3 / 2 pf exp(— p 2 /m p k^T). Averaging of Eq. (|A5j) with this distribution gives 
the cross section in the form (14T>1) with z/ pp given by Eq. (|47| , where A pp is calculated by averaging the p^-dependent 
factors of Eqs. (|A"4| and (fM|l . 



B. PHOTOABSORPTION DUE TO ELECTRON-PROTON COLLISIONS 

In the case of photoabsorption due to the electron-proton collisions, the initial and final states in Eq. (jAlj) are the 
continuum states of the hydrogen atom described by Eq. ([6|) . The wave function of the relative electron- proton motion 
can be written as 



ip{r) = ipo(r) +ipi(r), ip (r) 



e ik z 

~7l 



(Bl) 



where t^o(r) describes free motion with z-component of the relative momentum hko, L is the normalization length, 
and V'i( T *) is a perturbation due to the Coulomb interaction. Let us apply one-dimensional Fourier transformation 



L/2 



L/2 



' -0( r ) dz. 



(B2) 



In the limit of L — ► oo, we will have ipo(r±, fc) — > (2n/L) S(k — ko) Q ns (r±). Let us use expansion (TT2]) for ipi(r); then 

Mr ± ,k) =J2^'s'{k)^ ns (r±). (B3) 



This is equivalent to replacing g n > s '(z) by L x / 2 exp(ifcoz) 5 nn i5 ss i + g n i s >(z) in Eqs. (|12[) and (|15p . Then, applying the 
Fourier transformation to Eq. (|15p with t*o = r c and treating ipi as small perturbation, in the first approximation we 
obtain 



[(h 2 /2^) (fc 2 - fc 2 ) + E^, s , - E^ s ] ~g n , s ,(k) = -L- 1 V ns , n , s .(r c , k - fc ), 



where 



Vns,n' s' {j'c-} fe) 



Using Eqs. (A3)— (A10) of lPotekhinl <|1994D . we can convert V n 
where 



Vns.n' s'(r c ,z) dz. 

s,n's'{r c , k) into 
, v2a m k 



V2 



n p+ n P min(n p ,i) 

W' S '(p^) = (-i) (|shs+|s ' | - s ' )/2 E E °- 

/— m=max(0,2 — n' ) 



and, assuming q > 0, 

s „i 
W0, S ,0,s+g(/O, X) = V ,s+q,0,s{P, X) = P q E a msqP 2m / exp 

m=0 70 



x 2 1 -i 



p 2 i 



' dt. 



(B4) 
(B5) 

(B6) 
(B7) 

(B8) 
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In Eqs. (JB7J) and (TB8i we have denned 



n+(s-\s\)/2, n'=n' + (s'-\s'\)/2, 



s ~=(| s | +s + | s '|- s ')/2 + Z, 
g / = (M-a + l a / |+0/2 + ^ 
y n! (n + s)\ 



(B9a) 

(B9b) 
(B9c) 

(B9d) 



ml (n — m) \ (m + s)\ 

[rip and n' p are the radial quantum numbers of the Landau functions — e.g. jLandau &: Lifshitdll976f ). 

For fixed quantum numbers rif,Sf, and a fixed sign of the z- projection of the relative momentum kf of the final 
state, dvf in Eq. (|A1|) equals Ldkf/2n = (L/2n) (fj,/h 2 \kf\) dEf. Therefore, the cross section of photoabsorption 
by an electron- proton pair with initial quantum numbers rii and Si, longitudinal wave vector fc^, and transverse 
pseudomomcntum K± is 



a(ki, K±,rii, Si,uj) 



E 



^ |e.(/b cff |*}| 2 



nf,Sf ,signfc f 



n?\k f \ 



Here, the sum is performed over those n f and s f which are permitted by the energy conservation law, 



Ei 



h 2 k 2 f 

} - = 

2/i 



h 2 kf 

2/i 



hio. 



(BIO) 



(Bll) 



A general expression for j D g has been derived by Potckhin fc Pavlovl (|l997). In the dipole approximation, it reduces 
to 

tv n 



Jeff = e 



m m r 



where 7r is defined by Eq. ©, and 



II = p -B x r. 

y 2c 



(B12) 



(B13) 



The circular components of operators tt and II, n±i — (n x ± YEy)j\f2 and II±i = (LTj, ± illy) / v2, transform one 
Landau state |n, s)±, characterized by the function & ns (r±), into another Landau state, 4 

ih 



7T+l|n, s)± - 

K-i\n,s)± 
II + i|n, s}± - 



ift 



V« + l |n + 1, s - 1)_ 
|n — 1, s + 



ih 



Vn + s \n, s — 



II_i |n, s)j_ = — V« + s + 1 |n, s + 1) 



(B14a) 
(B14b) 
(B14c) 
(B14d) 



Since e • j = e+ij'—i + eoio + e -i i+ij the matrix element with 7r + i and n + i contributes to <j_i, and vice versa. 
In the first Born approximation, using Eqs. (|B4|) . (|B11|) . (|B12p .and (|B14|) . we obtain: 



(f\jes\i)o 



L fiuj 



(ki-kf)Vr 



TLiSiU f S f 1 



(B15a) 



(/b'eff|i)-i = - 



yjnf + lV niSi ,nf+l,Sf- 



TliV ni - 



l,Sj+l, nfSf 



m c (lo + lu cc ) 

y/llf + Sf V ni Si,nf,Sf-l — V n i + s i + 1 Kii,st+l,n/s/ 



+ 



(/lieffK) 



(B15b) 



?V V niSi ,n f -l,sf+l — V n i + 1 Vni+1,« 4 — l,n/s/ 



m c (w - W cc ) 



\/ n / + s / + 1 KiSi,n/,s/ + l - V n i + s i Vnusir 

m p (u + Lu cp ) 



(B15c) 



4 T he square-root factors in Eqs. (IB14ct and JB14dl ) were interchanged by mistake in corresponding Eq. (A3b) of Potckhin & Pavlovl 
Jl997h . 
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Here, for bre vity, y nsn > s > = V nsn's' (r c , kf - h). 

Equations (1B10|I and (1B15|) provide the partial cross sections for one electron-proton pair in a given state. Provided 
there are n e electrons per unit volume, the number of electrons interacting with a given proton and having ki in the 
interval dh and r c in the surface element d 2 r c is 



Here, 



dNi = n c L d 2 r c F\\ (h) dh T n 



T\\{ki) = h{2irnk B T)- 1/2 exp{-h 2 k 2 /2p,k B T) 



is the Maxwell distribution of the continuum states over ki , and 

Fn iSi = (1 - e~ 0e ) (1 - e~ Pp ) exp [-n^ e - (n* + Sl )(3 p ] 
is the Boltzmann distribution over rii > and Sj > — rij. Thus the total cross section is 

/OO POO 
^(ki)dki / KxdK±a(ki,K±,ni,Si,uj). 
Jo 

For every polarization, let us write er(cj, T, B) in the form of Eq. (|4"9"| . Then 

Q /■ OO 1 

A « = iE^ E / ^exp[-(/3 c + /? p )^ 2 /4 ] <sSin/s/ (^_^), 

n/,s/,signjnr/ " I /I 

where, taking into account Eq. (|B11[) . we have 

h 2 = *? + 4(n, - n/) + 4—( Sl - s f ) + 4^ ^, 

and the functions w" iS . nfSf (xf — in the integrand are defined according to Eqs. (IB6j) . (IB15|) . and (|B19|) as 

f°° 2 

Jo 

2 



+ 



(1 + m c /m p ) 2 J 
m c uj — uj, 



pdp 



a' ''„.-. „'-l,s'+l(/>, x) - V / n+T V„+l, s -l,n's'(P) x) 

2 



[Vn' + s' + lv nsn i tS ' + x{p, x) - Vn + sv n ,s-i,n'8'(p, 



(1 + m c /m p ) 2 J Q 



pdp 



Vn' + 1 Vns,n'+l,s'-x{p, x) - VnV n -l,s+l,n's'(p, 



m e co + uj, 



ii lq uj ^ ujcc r / — 7 / \ i / 

H \yn' + s' v nsn ^ s '- 1 (p, x) - vn+ s + li>n,s+i,n's'(p, 



m D a; — uj, 



cp 



(B16) 
(B17) 

(B18) 
(B19) 

(B20) 
(B21) 
(B22a) 

(B22b) 

(B22c) 
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